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Abstract. Selected highlights of the theoretical developments reported at the 2004 
Quark Matter conference are discussed, with emphasis on open issues. 



1. Introduction 

High-energy lieavy-ion collisions (HEHIC's) aim at creating the Quark-Gluon Plasma 
(QGP), together with a thorough investigation of the phase diagram of strongly 
interacting matter. The 17th International Conference on Ultrarelativistic Nucleus- 
Nucleus Collisions constituted another milestone towards these ambitious goals. 

From the theory side, the objective is to develop a coherent and consistent 
description of the relevant processes transforming the incoming nuclei into the finally 
measured (hadronic and electromagnetic) particle numbers, spectra, and correlations. 
The combination of different approaches is inevitable; e.g., the evolution of bulk matter 
might be cast as a succession of color-glass initial state, hydrodynamic expansion of 
ultradense phases, and transport simulations of moderately dense matter. To interface 
these building blocks, knowledge on the relevant microscopic degrees of freedom, their 
interactions and spectral properties is mandatory. The latter also hold the key to relating 
observables to (pseudo-) order parameters of the phase transition(s) to be identified. A 
further challenge is posed by the transition regime from bulk matter to perturbative 
spectra, as emphasized in Wiedemann's theoretical overview talk P]. 

Along these lines, I review below progress that has been reported on theoretical 
approaches to the QCD phase diagram (Sec. |21), on properties of bulk matter as deduced 
from RHIC runs 2+3 (Sec. and on microscopic probes expected to be scrutinized in 
run 4 and the recent SPS run (Sec. EJ. Final remarks are made in Sec. El 



2. QCD Theory at Finite Density and Temperature 
2.1. Color superconductivity 

The richness of the phase diagram at high densities and small temperatures continues 
to provide for new phases, especially under conditions relevant for compact stars such 
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as color and charge neutrality. The key quantity at moderate quark chemical potential, 
400 MeV, is the (constituent) strange quark mass, m^*\ Rajagopal elucidated that, 
reducing /ig from large values, the first transition is probably from the u-d-s symmetric 
CFL ground state into a "gapless" CFL phase j2]: the asymmetry induced by 
trifurcates the gap value and induces "blocking" regions on the d- and s- Fermi surfaces 
where their pairing is inhibited. One of the open questions is whether a 2SC state {u-d 
pairing only) exists in nature, which essentially hinges on the condition fig < m* being 
met in the quark phase, thus depending on nonperturbative forces responsible for chiral 
symmetry breaking (xSB) in the vacuum. 

Applications to neutron star phenomenology were discussed by Reddy, who argued 
that star radii of 8 km or less clearly favor a quark- matter equation of state (EoS). 
Much like in HEHIC's, electroweak emission is a promising signal of in-medium effects, 
probing the low-lying excitations (Goldstone bosons, (di)quark quasiparticles) of the 
matter in, e.g., supernovae explosions or long-term cooling 

2.2. Chiral vs. Deconfinement Transition 

A longstanding problem associated with the QCD phase diagram is the relation between 
chiral and deconfinement transitions. Current lattice calculations indicate that in real 
QCD the small but finite current quark masses, mg, lead to a finite-T transition which 
is not of first order. However, the peaks in the susceptibilities pertinent to order para- 
meters at rriq^oo and mq=0 - Polyakov loop (L) and quark condensate (qq), respectively 
- precisely coincide jBJ. Fukushima suggested [S] this feature to arise through a mixing 
of the scalar "soft" modes that become massless at the second order endpoints in 
the T-rriq plane, cf. left panel of Fig. ^ For the chiral transition this is the usual 
"cr" -meson, whereas for deconfinement he argued it to be the electric glueball field. 
Their mixing ensures the coincidence of the (pseudo-) transitions also in the cross-over 
region. It is not quite clear, however, why this mechanism is not operative for quarks in 
the adjoint representation where lattice QCD (LQCD) observes two distinct transition 
temperatures Similar issues were addressed by Mocsy [TUj . 

2.3. Lattice QCD Results 

One can close the gap between the two endpoints C and D in Fig. [T] (left panel) by 
moving to finite fig. Recent lattice computations have pursued this by an expansion in 
fiq/T, as discussed by Redlich and Karsch [HI [T2]. Fig. [T] (right panel) shows results 
for the quark number susceptibility, which indeed develops a rather sharp maximum 
(signaling a rapid change in baryon density) at yUB=3yUg— 500 MeV, quite reminiscent to 
independent determinations of the critical endpoint [T2j ■ This puts the latter into the 
realm of the future GSI facility, the physics of which was reviewed by Friman 

Further progress on finite-T (quenched) LQCD results for charmonium properties 
was reported by Karsch, Asakawa and Datta. From the color-singlet free energy, 
Fi = Vi — TS, the pertinent cc potential has been extracted which appears to support 
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Figure 1. Left panel: first order lines (solid) and critical endpoints for chiral (lower, C) 
and deconfinenient (upper, D) transitions in the T-rriq plane for 2+1 flavor QCD [71|S]. 
Right panel: quark number susceptibilities for different chemical potentials 

bound states for a range of temperatures above This corroborates the analyses of 

spectral functions exhibiting resonance/bound states of low-lying charmonia {rjc, J/ip) 
surviving up to ~2Tc fl(i|ll7j. which could have important consequences for charmonium 
production at RHIC, see Sec. 14.31 below. 

3. QCD Bulk Properties at RHIC 

In this section, I will essentially follow a chronological order in discussing the properties 
of matter as evolving in a HEHIC, i.e., from the initial state (Sec. 13.111 via the dense 
phases through the phase boundary (Sees. 13.21 13. 3|) until thermal freezeout (Sec. 13. 4|) . 

3.1. Color-Glass Condensate (CCC) 

An important experimental finding since the previous Quark Matter conference is the 
absence of a suppression of (moderately) high transverse-momentum {pt) hadrons in 
200 GeV d-Au collisions at midrapidity (77=0), which, however, turns into a suppression 
at forward rj ^H|; cf. left panel of Fig. |21 Whereas the former confirmed strong 
suppression of high-pj hadrons in Au-Au as an energy loss effect in the produced 
hot/dense medium (see next section), the latter might have provided first explicit 
evidence of gluon saturation in the low-x part of the Au nucleus' wave function, as 
advocated by Jalilian-Marian, Venugopalan jTHj and Kovchegov (see also Ref. [201 )• In 
addition, the increase of suppression with centrality of the d-Au collision is qualitatively 
in line with the CGC j25- Definite conclusions have to await quantitative comparison 
with data incorporating also more mundane mechanisms, as stressed by Accardi |22j . 
E.g., the HIJING event generator |2S1 satisfactorily describes the rapidity spectra in 
d-Au [25, as well as the h'^ /h~ asymmetry in the forward region, indicative for valence- 
quark fragmentation. The former implies that the low-pj regime is not an unambiguous 
signal of CGC. The critical region will thus be around the estimated saturation 
momentum of Qs— 2 GeV (for 7^=3 at RHIC). Forward charge asymmetries further 
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Figure 2. Left panel: rapidity dependence of RdAu{pt) in 200 GeV d-Au JHI- Right 
panel: STAR data |2H| and jet quenching calculations [2^1 for azimuthal correlations. 



complicate the interpretation of the nuclear modification factor RdAu which uses a p- 
p reference as denominator. Rapidity energy loss of projectile valence quarks can be 
expected to induce a decreasing centrality dependence of RdAu, too. 

One should note that if the CGC is indeed operative at moderate pt in forward d-Au, 
it will refiect itself also in bulk hadron production systematics in Au-Au at midrapidity. 
Using pertinent initial conditions for a hydrodynamic evolution leads to a satisfactory 
description of the observed dN/dy distributions, as shown in Nara's talk |25j . 

3.2. Opacity and Thermalization 

The main evidence for the production of thermalized matter at RHIC, at energy densities 
well above the critical one, resides on the quantitative success of parton energy-loss and 
hydrodynamic calculations. 

As discussed by Vitev and Barnafoldi, using well calibrated perturbative QCD 
probes in connection with mostly medium-induced (non-abelian) gluon radiation allows 
to infer energy densities e~20 GeV/fm^ in the early stages of 200 AGeV central Au-Au 
collisions. In addition to correctly predicting the factor of 4-5 reduction of the nuclear 
modification factor, RAA{pt), together with its flatness up to the highest currently 
measured pt—11 GeV, the calculations reproduce azimuthal correlations, especially for 
the away-side jet, such as its (slight) broadening in d-Au its gradual disappearance 
with centrality in Au-Au |22i (cf. right panel of Fig.|21), as well as its preferably in-plane 
reappearance in semi-central Au-Au (in the short direction of the "almond" ) . 

The relevance of "hadronic" quenching was addressed by Greiner. Whereas Lorentz- 
dilation precludes noticeable impact on formed hadrons, absorption of colorless "pre- 
hadrons" (emerging from hard-scattered valence quarks color-neutralizing on timescales 
t^Ejet/Q'^^l/pt |2n|) on surrounding (pre-) hadrons may be signiflcant. However, the 
implementation of this approach in a transport model [SOI tends to overestimate high-p^ 
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Figure 3. Left panel: spectral functions for light hadrons as extracted from quenched 
LQCD calculations [SJ- Right panel: schematic dependence of light hadrons masses 
across the phase boundary, with S- and P-wave dissolution points in the QGP |33| . 



suppression at SPS, while underestimating elliptic flow at RHIC. 

Ideal hydrodynamics is remarkably successful in describing the collective evolution 
of ~ 99% of all produced particles at RHIC as illustrated in Hirano's talk (corrections 
away from midrapidity and due to viscosity were discussed by Heinz and Teaney, 
respectively). The inherent early thermalization time, ro~0.5 fm/c, required especially 
for second and fourth order azimuthal asymmetries, V2{pt) and V4^{pt), however, cannot 
be accounted for in terms of perturbative partonic rescattering cross sections, which 
are much too small. Here, progress could lie in the realization that for moderate 
temperatures (1-3 Tc) the QGP may support q-q and/or g-g bound states (resonances), 
as evidenced by recent LQCD calculations of "hadronic" spectral functions jHI], cf. left 
panel of Fig. El Shuryak and Zahed ^32j suggested the color-Coulomb interaction 
as the underlying force, with the running of ag only limited at the screening scale 
rriDebye- Browu argued that, close to Tc, additional instanton-molecule (/M) induced 
interactions (augmented by RPA resummations and spatially small Coulomb wave 
functions) become important to bind thermal quasi-quarks into (almost) massless 
pions being continuously connected to the hadronic phase (even though standard 
estimates of IM coupling strengths are rather small jSl]), cf. right panel of Fig. El 

Resonant rescattering of partons is a promising mechanism to explain the short 
thermalization times required by hydrodynamics, provided the resonance correlations 
themselves build up sufficiently fast. 

3.3. Hydrodynamics vs. Quark Coalescence 

Relativistic hydrodynamics, employing an equation of state consistent with predictions 
from LQCD, provides an excellent description of low-p^ hadron spectra. Implementing 
jets and their energy loss into a hydrodynamic evolution, Hirano and Nara confirmed 
that the transition momentum from collective to perturbative production increases with 
mass (2, 2.6 and 3.5 GeV for tt, K and p in central Au-Au, respectively) |35j. However, 
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even the hydro protons appear to "run out of steam" too early, in that the p/n ratio 
in the hydro+jet model underpredicts the experimental data above pt— 3 GeV 
(similarly, the experimental A/K ratio recovers the perturbative value only close to 
Pi~6 GeV flZI)- Furthermore, the elliptic flow of all thus far identifled hadrons appears 
to follow a constituent-quark scaling, i.e., a universal "partonic" V2{pt/n)/n with n=2 
(3) for mesons (baryons). Both the universal V2 and the enhanced baryon-to-meson 
ratios at intermediate pt are readily explained within a quark coalescence picture at 
hadronization, as shown in the talks by Fries, Hwa and Molnar. Whereas the model 
of Fries et al. jHHI focuses on coalescence of thermal partons, Greco et al. [2^] and 
Hwa et al. jlO] also allow for recombination with minijet partons. Whether the latter 
will be able to explain the observed nearside correlations in the 4 GeV regime (recall 
right panel of Fig. |2)) - which appear unmodifled from p-p to central Au-Au - remains 
an open question at present. In fact, in the approach of Hwa et al. PU^, also in p- 
p collisions a signiflcant contribution is assigned to coalescence. As for the distinction 
from hydrodynamic models, the elliptic flow of mesons (with a mass similar to baryons 
but only n=2 constituent quarks) is an ideal observable. First data on the ratio of central 
to peripheral prspectra, R^pi^pt), seem to follow meson systematics.|| 

In view of the LQCD spectral functions (Fig. E)), a logical extension of coalescence 
models could include hadron formation even above Tc, e.g., by solving suitable rate equa- 
tions, reminiscent to recent calculations in the charmonium context ^] (see Sec. 14. 3j] . 

3.4. HBT "Puzzle" 

New suggestions have been presented at this meeting for the problem of the HBT data, 
especially the smallness of the "out"-radius, R^^^ = D{xout,Xout) — '^D{xout, Ptt) + 
D{(3tt,pt) {D: variance). Kapusta jl2] and Wong argued that keeping track of 
quantum phases in the hadronic rescattering could preserve memory on the initial source 
size, thereby also addressing the approximate constancy of the radii with collision energy. 

We recall that the AMPT transport model ^1] can account for the measured radii, 
with a positive Xout-t correlation, (partially) related to decays of long-lived resonances 
(especially uj's), as a key ingredient to reduce Rout- The question remains where the 
discrepancy to other transport models ^| resides, which should be rather similar in the 
treatment of the late stages, thus pointing to differences in earlier stages. 

4. Microscopic Probes of QCD Matter 

With the bulk properties of the produced matter at RHIC being reasonably well assessed, 
the next challenge is to determine its microscopic properties. In this section, I will focus 
on 3 complexes of observables that are expected to serve this goal, proceeding from the 
relatively dilute freezeout to the hot and dense phases. 

II One should note that hydrodynamics ceases to be appHcable m peripheral collisions. 
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Figure 4. Left panel: model fit |i48| to STAR ,47^ tt+tt invariant-mass spectra with 
a p-mass shift of -50 MeV. Right panel: direct photon sources at RHIC. 



4.1. Resonance Spectroscopy 

Short-lived resonances are a promising tool to deduce in-medium modifications of their 
spectral shape close to thermal freezeout through invariant-mass spectra of their decay 
products, e.g., p° ^ tt+tt- or A++ ^ p7r+ For (/z^, /i^, T)jo=(370,90,110) MeV, 
the pion (anti-/baryon) density is, in fact, still appreciable, ^^=0.65^o (Qb+b^^-'^^Qo): 
with ^)o=0.16 fm~^. The left panel of Fig.|3]shows a model fit to STAR data [Uj from 40- 
80% central Au-Au collisions employing empirical vrvr phase shifts in a sudden breakup 
scenario as discussed by Florkowski |1H]. He inferred a p-mass drop of about 50 MeV, 
similar to previous findings [101 EOl- Broadening effects, Bose-Einstein correlations 
and nonresonant background contributions may reduce this value significantly |5H I52j . 
Also, a comprehensive understanding of appreciable mass shifts as extracted from high- 
multiplicity p-p or e'^e~ reactions is pending. Preliminary data for the A"'"^ resonance 
show a slight increase of its mass, together with a significant broadening, when going 
from peripheral to central 200 AGeV Au-Au collisions |46j . 

Concerning the resonance yields, measured p/tt, Jq/tt and A/p ratios ^| are at 
least a factor of 2 larger than equilibrium values at thermal freezeout j3T] , which could be 
suggestive for multiple emission necessitating a rate equation approach [IHIEDIEII, e.g., 
2-3 generations of emitted p, fo and A's. Such a scenario can simultaneously account 
for the observed K* (892) / K ratio [46 - which is reproduced by the equilibrium value 
at Tfo WD^ - since ri^.~2-3 Tpj^^i^. 

Towards the phase boundary, chiral symmetry restoration eventually dictates sub- 
stantial medium modifications in terms of degenerate spectral functions for chiral part- 
ners, e.g., TT-o", p-ai, A^-A^*(1535). The dilepton enhancement found at SPS [Sni (to be 
scrutinized by NA60 [33] , as well as PHENIX at RHIC) may well be a related signal in 
the p—>-ee channel [SS] ■ On the one hand, this poses the challenge to be reconciled with 
chemical-freezeout models jSHl EZ| (usually based on vacuum masses), e.g., using self- 
consistent thermodynamic potentials [SHI- the other hand, resonance spectroscopy 
may offer direct access to pertinent chiral partners, e.g., via ai^ir'y or N*^r]N. 
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4-2. Electromagnetic (EM) Radiation 

The penetrating character of EM emission makes it a versatile probe of all collision 
stages. It should be emphasized that photons and dileptons are (kinematic) facets of 
the same object, i.e., the electromagnetic current correlation function (for M^=0 and 
Mee>0, respectively). Gale in his talk advocated that both can be used to tag initial 
energies of jets jSHl IHOl; the latter also radiate electromagnetically in their passage 
through matter, which, in particular for 7's, could constitute a major source for pt—2- 
6 GeV ini], possibly outshining both pQCD and thermal radiation. A characteristic 
signature of this contribution is its negative elliptic flow. EM radiation off jets is part of 
what is more generally denoted as "pre-equilibrium" sources, which have recently been 
evaluated in a parton cascade simulation as presented in Bass' talk jHll- The pertinent 
photon spectra for RHIC are confronted with hydrodynamic jHSj as well as thermal 
fireball calculations [Hll in Fig. El (right panel), also suggesting that pre-equilibrium 
yields dominate above pt—2 GeV. The largest component in the parton cascade is due to 
rescattering of secondary partons (quantum interference (LPM) effects are not included). 
The two thermal calculations agree rather well, with hadron-gas radiation prevalent up 
to pt—1 GeV, where QGP radiation takes over. It is also noteworthy that chemical off- 
equilibrium in a QGP affects the thermal photon yields very little, i.e., undersaturated 
parton distributions are essentially compensated by higher temperatures [HH] . 

WA98 at SPS reported new data on low-p^ (100-300 MeV) direct photons ^3] 
which exhibit a pronounced excess of a factor of 3 or more over thermal fireball calcu- 
lations ini] . The latter employed a recently improved assessment of hadronic rates which 
simultaneously describe CERES low-mass dilepton data. If confirmed, the \ow-pt photon 
enhancement could point to appreciable in-medium effects in the late hadronic stages 
currently not accounted for, e.g., Bremsstrahlung from tttt scattering via in-medium 
softened "a" -mesons (a possibly related softening might have been observed in low- 
energy vr-induced vr- production off nuclei, hA^tiitA, see, e.g., Ref. pB]). One should 
note, however, that the CERES data at very low mass (Mee<100 MeV) do not leave 
much room for strong additional sources |53j . 

4-3. Charm (onium) 

The intermediate scale of the charm-quark mass renders charmed hadrons another 
valuable messenger of QGP properties. The first step is to establish a baseline for 
their primordial spectra. A CGC calculation presented by Tuchin predicts a 40% 
decrease in dN/drj for open-charm when going from ?7=0 to 2 in central Au-Au, which is 
essentially a consequence of the saturation scale changing from Qs<mc at ri=0 (implying 
A^co«-scaling) to Qs>mc at ri=2 (entailing A^^pa^r scaling). On the other hand, in the 
color-dipole approach of Raufeisen's talk [HHI, a 25% shadowing is found for midrapidity 
yields which does not significantly vary for r]<2. 

Recent PHENIX data for "non-photonic" single-electron pt-spectra jnE] (ascribed 
to open heavy-flavor decays) in central Au-Au are consistent with the "null-effect" on 
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Figure 5. Left panel: PHENIX elliptic flow data for "charm-like" single electrons |7T] 
compared to coalescence model predictions for D-^evX (72,. Right panel: PHENIX 
J/^ data IZHl compared to a kinetic rate equation approach |^ . 



the c-quark energy loss, see Djordjevic's talk [70]. However, first PHENIX data on the 
elliptic flow of single-e^ jTT] are suggestive for a non-vanishing signal, cf. left panel of 
Fig. El The plot also includes coalescence model predictions [72] for parent D-mesons 
formed by recombining thermal light quarks (including (elliptic) flow as determined by 
pion data [HI]) with c-quarks either from PYTHIA (representing no rescattering, dashed 
line) or being thermalized in, and flowing with, the bulk matter (solid line). No deflnite 
conclusions can be drawn yet. Note that, although the two opposite scenarios of no 
c-quark reinteraction vs. thermalization lead to very similar single-e^ pt-spectra up to 
~2 GeV [ZS], pertinent c?z-lepton invariant-mass spectra should be quite different, since 
the back-to-back character of hard production (implying large M^e) becomes randomized 
in the thermal case (reducing the average Mee). 

Thermalization of c-quarks has important impact on charmonium production. With 
10-20 cc pairs in central 200 AGeV Au-Au, and the realization from LQCD that J/ip 
states possibly persist up to 2Tc, their regeneration in the QGP becomes plausible, e.g., 
via the inverse of gluon dissociation reactions, J/ip+g ^ c+c+X^. In particular, 
thermalized c-quarks enable the description of these processes in a kinetic theory 
framework via simplifled rate equations of type dN^/dt = V^{N!^ — N^), where the 
reaction rate represents the width of in-medium J/ ip spectral functions, and thus, in 
principle, is directly amenable to LQCD calculations. In addition, if, as to be expected, 
the open-charm number is conserved in the course of an A- A reaction (i.e., determined by 
primordial N-N collisions), the equilibrium level of charmonia is sensitive to open-charm 
masses [H], which can also be determined in LQCD. Grandchamp in his talk showed 
solutions to the kinetic equation including essential features of LQCD [H], cf. right 
panel of Fig. El Most of the J/ip^s in central Au-Au at RHIC are indeed regenerated in 
the QGP, with the band indicating sensitivities to in-medium open-charm masses. The 
observation of a signiflcant number of J/^/^'s at RHIC would be another step towards 
ascertaining the notion of hadronic resonance states in the (nonperturbative) QGP. 

^ This may be viewed as a generalization of statistical production at Tc [T^, see also Refs. [7511751177] 
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5. Final Remarks 



The search for the QGP has reached a critical phase: the combination of current RHIC 
data with theoretical analyses shows that the matter created in central Au-Au collisions 
is (i) very dense, (ii) thermalized, and (iii) probably nonperturbative. Further combining 
(i) + (ii) with the lattice QCD value for the critical energy density, one is lead to conclude 
that the QGP has indeed been produced at RHIC. However, to claim such a discovery, 
one would like to have a deeper understanding of the nature of that phase, as encoded 
in (iii), i.e., what the relevant (microscopic) degrees of freedom and their interactions 
are. In this respect, upcoming charm(onium) and electromagnetic probes (photons, 
dileptons) can be expected to provide valuable insights. 
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